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In this work we calculate the low-energy eflFective action for gravity with torsion, obtained after the 
integration of scalar and fermionic matter fields, using the local momentum representation based on 
the Riemann normal coordinates expansion. By considering this expansion around different space- 
time points, we also compute the nonlocal terms together with the more usual divergent ones. Finally 
we discuss the applicability of our results to the calculation of particle production probabilities. 
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^ I. INTRODUCTION 

o . 

' The absence of a quantum theory of gravitation is probably one of the most important open problems in theoretical 
physics. In fact all of our knowledge about gravitation refers only to its classical aspects virhich arc well described by 
^sj ' General Relativity and other generalizations of this theory with the same low-energy limit. On the other hand, the 
other interactions (strong and electroweak) are well described, even at the quantum level, by the Standard Model, 

■ as it has been confirmed in detail in the last years at the Large Electron-Positron Collider (LEP) and many other 
experiments. Moreover it seems reasonable to think that there may exist some physical regime where the gravitational 
field could be treated classically whereas matter is quantized. In fact this is the only regime that can be studied by 
now without making extra hypothesis beyond what has already been checked experimentally. In this very conservative 
approach one considers, on the one hand, the quantum gauge and matter fields propagating on a curved space-time 

T— I ' and, on the other hand, the dynamics of the classical degrees of freedom associated to the space-time. This is the 
t — ■ ', so-called semi-classical approximation. It is clear that space-time curvature or torsion (see for a review) can affect 
0^ ■ the quantum fields in different ways. In addition, space-time is also affected by the presence of the quantum matter 
' fields, as it happens, for instance, in the electromagnetic field case, where fermion loops modify the electromagnetic 
. dynamics by means of the vacuum polarization and other effects. 

In order to study these new effects due to matter fields loops, it is specially appropriate the use of the gravitational 
D . effective action (EA). This EA is obtained after integrating out the matter fields. In general the EA will be a nonlocal 

■ and non-analytical functional of the metric and the connection. The exact expression for the EA, obtained after 
integrating out the matter fields W^[fl';iy] is not known for arbitrary space-time geometries. However there are several 
techniques that have been proposed for its approximate calculation, namely: perturbation theory, in which the metric 
tensor is divided in a fiat metric and a small fluctuation. The main difficulty of this method is that general covariance is 

^ . explicitly lost ^-|^ . The Schwinger-DeWitt proper time representation allows us to obtain a covariant asymptotic 
■ ■ ' expansion for the EA, but as far as it is local, it does not allow to describe particle creation processes. These difficulties 
are intended to be solved by means of the partial resummation of the Schwinger-DeWitt series Q . There are also 
perturbative methods that respect general covariance, such as the so called covariant perturbation theory Finally 
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the local momentum representation 
the advantage of combining the usua 



|10j] is based on the Riemann normal coordinate expansion. This technique has 
flat space-time methods with formally covariant expressions. It has allowed to 
calculate the divergent local parts of the one loop effective action for scalar and fermionic theories. The main aim of 
the present work is to find a representation of the nonlocal finite parts of the effective action in this formalism. 

Once we know the EA, we have all the information concerning the semi-classical gravitational evolution. The 
corresponding equations of motion will modify the Einstein field equations taking into account quantum effects. 
Moreover, the EA could have a non-vanishing imaginary part, which can be interpreted as the particle production 
probability §,0. 

In this work we show our results concerning the computation of the one-loop EA after integrating out scalar and 
fermionic fields by using the local momentum representation. Our computation includes not only the divergences but 
also the nonlocal finite terms that can lead to instabilities of the classical solutions by particle emission. As we will 
show, the use of normal coordinates allow us to obtain a useful representation of the nonlocal form factors. Such 
representation has been succesfully applied in a recent work to obtain particle production in different cosmological 



contexts in a remarkably easy way |11|. In addition, the results will shed some light about the boundary conditions 
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on the metric tensor which permit a definition of the effective action. 

The work is organized as follows. In Section 2 we do a brief review of the Euler-Heisenberg effective lagrangian for 
Quantum Electrodynamics (QED). This model will be a guide for the calculation of the corresponding gravitational 
EA. In Section 3 we present the method to generate the derivative expansions of the EA by means of Riemann normal 
coordinates. Applying such method to the scalar theory in the presence of gravitation, we obtain the divergences as 
well as the finite nonlocal pieces of the EA up to quadratic terms in the curvature. In Section 4 we do the same with 
fermionic matter fields, obtaining the corresponding gravitational EA up to quadratic terms in the curvature and in 
this case also in the torsion. In particular we apply our results to the Standard Model particle content. Finally Section 
5 contains the main conclusions of this work an a brief discussion on their possible applications to the computation 
of particle production probabilities. We have also included an Appendix containing the dimensional regularization 
formulae and some normal coordinates expansions used in the text. 



II. THE EULER-HEISENBERG LAGRANGIAN 



The historical origin of the semi-classical EA can be traced back to the Euler-Heisenberg lagrangian for QED . 
When the momentum p of photons is much smaller than the electron mass M, the one-loop effects, such as vacuum 
polarization, can be taken into account by adding local non-linear terms to the classical electromagnetic lagrangian. 
Consider the QED EA given by: 



(1) 



where as usual p ~ ^^{d^ — ieA^). Using dimensional regularization (see Appendix), it is possible to find the following 
expression up to quadratic terms in the photon field: 

W[A] =-\j d'^xF^.F^'^ - zTrlog((z p~M + ze)) = -\ J d'^^F^.F^'^ + i ^^Tr[(z^- M)-^ 



k=l 

2e2 „„„ /2M2 If 



-F^'^F.^ — -Ff'^Fa^ + - — -F^"" + - 1 - 2 F(-n; M^) F^u 

4 '"^ 3(47r)2 (Att)^ ^3 □ 6 V □ ' ^ ^ 



-0(^4) (2) 



where A = — \og{M'^ / ^'^) , with = 2/e + log47r — 7 and 7 ~ 0.577 is the Euler constant. We have performed 
the formal Taylor expansion of the logarithm and used the expression: 

F{-n;M')F,^{x)= J d'y^e^'^(^-y^F{p';M')F,^{y) (3) 

with: 

/•I / „2 \ 

(4) 



In a similar way, the inverse operator l/D can be defined with the usual boundary conditions on the fields as: 

-^F,A^) = / d4y^e*P(-^)^-F^,(y) (5) 
— □ J (27rJ^ p'' + le 

The expression for the EA is nonlocal and has a regular massless limit. In fact, for small p compared with M, the 
Mandelstam function F{p^; M^) behaves as: 

F{p';M') = - log ( f . ) + 0{M^) (6) 

From we can see that the only contributions in the massless limit are those coming, on one hand from the A 
factor and, on the other hand, from the Mandelstam function. Both logarithmic contributions equal, up to sign, so 
that they cancel each other and we obtain: 

W[A] = J d'x V^^F^. - -l-^F^^^nu)F,^ + 0{A') (7) 
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where we have used the following notation: 

r(n) = TV, - log (^^) (8) 

to be understood as in the previous cases through the corresponding Fourier transform, with the ie factor as shown 
in (§. 

The EA (|^) allows us to derive in an exact fashion the photon two-point one loop Green functions. This, in turn, 
ahows us to obtain for example the vacuum polarization. The EA can be expanded as a power series in p^/A/^, and 
also in A to obtain the well-known Euler-Heisenberg local lagrangian 



III. INTEGRATION OF MATTER FIELDS IN A GRAVITATIONAL BACKGROUND 

Consider a scalar field in a curved space-time. The corresponding classical action is given by: 

= -\j d^x^cj, (□ + + ei?) (9) 

where □</() = g'^'^V ^idv4> = g~^^^d^ {g'^'^ \/9diy(j}) ■ The non-minimal term ^i? is included so that for to = and ^ = 1/6, 
the classical lagrangian is invariant under local conformal transformations. 

The EA for the gravitational fields that arises after integrating out real scalar matter fields is given by the following 
expression in Lorentzian signature: 



where Oa;j,(TO^) — {—Oy — — ^R{y) + ie)S^{x,y) with S^(x,y) being the covariant delta 5" (x, y) ~g ^^^{x)6{x,y). 
Therefore we can write: 

W[g^,] = ^log det 0(m2) = ^Tr logO(TO2) (11) 

Since in this expression we have only integrated the scalars out, the gravitational field is treated classically. Ac- 
cordingly, this EA is analogous to the classical action but including the quantum effects due to the matter fields. 
In addition, ( |l0[ ) is the generating functional of the Green functions containing only scalar loops and gravitational 
external legs. 

Let us now consider a fermion field propagating in a curved space-time with torsion . The corresponding classical 
action will be given by: 

" ^ \ (12) 



S = / d^xy/gip{i p - M)%P = / d'^x^'i}^ 



where S*^ = Cpi/A/^r'"'^ is the torsion pseudotrace, are the Levi-Civita spin-connection components and Saj, are 



the Lorentz group generators. The fermionic EA is given by 



^w[eX,A] = j \d%ljd%p] exp [i J d^xy/^ tpiip -M + ie)V^ = det(i p - M + ie) (13) 

where e denotes the vierbein, F the full connection with torsion and A collectively denotes the possible gauge fields. 
Accordingly: 

M^[e,f,^] = -ilog det(i ^ - M + ie) = -iTrlog(i ^ - M + ie) (14) 

In the massless limit, M = 0, the classical fermionic lagrangian is also conformally invariant. 

Since the above model does not posses self-interactions, the one-loop calculation is exact. This does not mean that 
it is possible to explicitly calculate the EA for an arbitrary space-time geometry. In those cases in which there is a 
high degree of symmetry, as in maximally symmetric manifolds, or in the so called conformally trivial situations, i.e, 
conformally flat manifolds and conformally invariant theories, it is possible to find the explicit form of the modified 
Einstein equations coming from the EA |0| . 
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In order to consider more general geometries, we will use an approximation scheme similar to the one used for 
the Euler-Heisenberg lagrangian. It consists in treating the curvature as a small perturbation. When we integrate 
massive fields out, this is equivalent to consider that the Compton wavelength corresponding to the massive particle 
is much smaller than the characteristic length scale of the gravitational field. In this framework, the expression for 
the EA will be an expansion in metric tensor derivatives over the particle mass. Such expansion will be generated by 
the normal coordinates expansions of the 0{m?) and {ip — M) operators. In the massless case, or if we are interested 
in the high-energy regime, it is possible to obtain an alternative expansion in powers of the curvatures (Riemann and 
Ricci tensors and scalar curvature), generically denoted 7?,. 



A. Riemann normal coordinates 



Let Xg be the coordinates of P in a given coordinate system and consider the set of geodesies passing through P, 
that we will write as a;"(r). We will choose the r parameter in such a way that a;"(0) = x^. Each of these geodesies 
will be characterized by the tangent vector at P, 



(15) 



and each point A on each geodesic by certain value of the r parameter. The Riemann coordinates of A are defined 
as = ^"r [Q. In a neighborhood of P where any other point A can be joined to P by a unique geodesic (normal 
neighborhood), the correspondence between the a;" and y" coordinates is one to one. 

When torsion is present and it is completely antisymmetric, the geodesic equation agrees with that obtained using 
the Levi-Civita connection and therefore, some given Riemann coordinates respect to the Levi-Civita connection will 
also be Riemannian respect to the connection with torsion ||l^,|l^ . By means of a linear real homogeneous coordinate 
transformation, it is possible to write the metric tensor at yo in the Minkowski form jy^^. The new coordinates are 
also Riemannian and they are known as Riemann normal coordinates. 

Let us consider the components of some tensor field at y, that we will assume to be analytic functions in a 
neighborhood of the normal coordinates origin j/o- From their Taylor expansion around we can obtain an expression 
that is written as a series in curvatures and their covariant derivatives. In paricular for the metric tensor componentes 
we get p|: 



1 



1 



+ 0(9^) 



y'^y^y'y' 



(16) 



Here, 0(9^) denotes terms with 5 or more derivatives and the indices in those tensors evaluated in j/o are raised 
and lowered with the fiat metric r\p^v. In the Appendix, we have written the expansions corresponding to the metric 
determinant and other useful expressions. It is important to note that (16) is a superposition of two kind of expansions: 
on one hand, terms are organized by the number of metric derivatives, but on the other hand, those terms with a 
certain number of derivatives can be classified according to the number of curvature tensors they have. 

It will also be useful to introduce the following relation: 2cr(x,a;') = 2/a2/", where the biscalar 17(0;, a;') represents 
half of the geodesic distance between the x and x' points and denotes the normal coordinates of the x point with 
origin at x' . On the other hand, d^a{x,x') is a tangent vector at x to the geodesic joining x and x' , whose length 
equals the geodesic distance between these two points and it is oriented in the — *■ x direction. In turn, d^'a{x,x') 
is tangent to the same geodesic at x' , with the same modulus and oriented in the opposite sense. 

In normal coordinates with origin at x' , according to the previous expressions, we can write: 



aaix,x') = -^a{x,x') = ya 



(17) 



i.e, ya are components of a vector tangent at the origin. 

The use of normal coordinates, apart from being basic to obtain the derivative expansions of the EA, allows us to 
work in momentum space in a similar way to the fiat space-time. Let us consider some scalar function f{x,yQ) with 
normal coordinates with origin at yo- We can define its covariant Fourier transform through ||l0|| : 



/(a;,yo) = 



/(fc,2/o)e" 



(18) 
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In a similar way we can introduce the covariant Dirac delta: 




As far as one of the delta arguments is the origin of coordinates, we will have S^{x, yo) = S{x, yQ)g ^^'^(yo) = Vo)- 
In the general case, (arbitrary arguments) there is also a covariant definition whose expression in arbitrary coordinates 
is given by 

S''{x,y) = g'/\x)g^/\y)A^/'{x,yo)A^/'{y,yo) J ^.g-^/^Cyo) exp(»fc^K(y, y„) - yo))) (20) 

where k^, <J^{y, yo) a-nd (Tf_i{x, yo) are tangent vectors at yo and A(x, x') is the Van Vlcck-Morcttc determinant, defined 
as: A(a;, x') = g~^/'^{x) det(— V^' Vi/(T(a;, x'))g^^^'^{x'). If we take yo as origin, we will have A(a;, yo) = g^^^^{x), this 
expression reduces then to: 

S'ix,y)^g-'/'ix) I -0^cM-^k,{x^-yn) (21) 
All these definitions are valid only in normal neighborhoods of the origin in which geodesies do not intersect. 



B. Derivative expansions 



In normal coordinates there is a privileged point around which we perform the expansions. In addition, the different 
curvature tensors are defined on the tangent plane corresponding to that point. This fact, together with the general 
coordinate invariance of the EA will allow us to obtain a covariant derivative expansion for the effective lagrangian 
around the origin. In the following we will discuss in detail the scalar case, although the procedure is the same for 
fields with different spin. 

Let us start with the scalar EA (11). Using the normal coordinates expansion for the metric tensor it is easy to 
split the operator Oxy{ni^) — {—^y — rn? — £,R{y) + ie)S'^{x, y) in a free part, that coincides with the flat space-time 
Klein-Gordon operator 



Axy{m^) = {-Ul-m^ + ie)5\x,y) 
with Dq = rj'^^d'j^dy, and the interaction part B^y that includes all the curvature dependence: 



(22) 



-lR\iyo)y''dl + lR''/piyo)yVdpy - myo) 

1 



20 



Ra;l3^iyo) 



— — R 

20 " ■^'^ 



(yo) 



20 



(yo) 



^RVp;^^yo) - ■^Rc.x{y^)RW{y^) + Y^^//3A(yo)i?\%(yo) + -^RcT pMR^iy^) + ■^Ro.p-,; iy^) 



45 



15 



45 



+ li?„^^%(yo)) y'^y^y-'dl - (^-^R^ ; ^.^^{y,) + l-R^^^^^[y,)R\'^ ^(^y,)^ y^y^y'y'^dldl 



i^R-Myo)y''y^ 



S"{x,y) + ... 



(23) 



We have only written the two and four derivatives contributions since terms with an odd number of metric derivatives 
are shown to be irrelevant for the final result. Therefore we have: Oxyijn?) = Axy(w?) + B^y 

We will also assume that space-time is asymptotically flat and this will allows us to discard total derivatives of the 
curvatures in the EA. We have included the covariant Dirac delta S'^{x, y) in the definition of the free operator so that 
we can use the covariant integration measure d'^xg^^'^{x). Taking all this into account we can write the EA as: 



W[g^.] = J 



= I <fxCeff{x) = -TrlogO(TO^ 



-Trlog(A-)-B) 



(24) 



As far as the calculations will be done in normal coordinates with respect to the yo point, any term in the effective 
lagrangian is evaluated in that point. For that reason the integration involved in the Tr symbol cannot be done 
immediately. This, in turn, keep us from expanding the logarithm, which only makes sense inside the trace. A way 
to avoid the problem consists in formally differentiating the EA with respect to w? |l^ : 
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In the last step the inverse oi A-\- B has been expanded. This can be done without taking the functional trace. 

The lowest order term in B is linear in the Riemann tensor and therefore it contains two metric derivatives. The 
term in (^5|) with two B factors will be 0{TI^) and will contain at least four metric derivatives, etc. Since we are only 
interested in the result up to 0(7i?), it will be enough to consider the first three terms in (p5|). 



C. Divergent parts 

In order to calculate the divergent parts of the EA, we consider the previous expansions. Since normal coordinates 
allows us to do momentum space integrations, we can use the flat space-time regularization techniques such as 
dimensional regularization. Our calculation method is based on the local momentum representation proposed in [ p^ , 
but adapted to the EA techniques. The results agree with the well-known expressions in p^ . 

First we find the scalar propagator A~y , given by: A^yAy^ — (5" (a;, z), where we have used the De Witt generalized 
summation convention for repeated indices. Using we can write this equation as: 

dW/'iy){-^^o - + / = (26) 

The G{k) function is easily obtained and finally the propagator reads: 

A-i ^ ^ r d*k ik(y-z) 1 (27) 

The first term in the EA expansion (pq) can be immediately evaluated and it reads: 



,-1 _ fs 1 ^_Til-D/i)^^ 

yoyo J fc2-m2 (47r)^/2 (m2)i--D/2 ^ 

with the notation dk — kfj,'^ /{2tt)^ . When i/q appears as a repeated subindex it must be understood that the 
integration in j/q has not been done. This final integration, that corresponds to the trace in (p5|), will be performed 
below in an explicit way. In the last step we have used the equation ( ^9[ ) from the Appendix. Performing the 
integration we have the lowest order term in the effective lagrangian: 

^2iyo) = ^m^(A+i) (29) 



647r2 V 2, 

where we have defined A = iV^ -t- \og{fi^/m^). As is well-known, this first term will give rise to the cosmological 
constant renormalization. The integration constant can be set to zero without loosing generality since it can be 
absorbed in the renormalization procedure as we will show below. 

Up to two metric derivatives, only the second term in ( p5| ) contributes: 

iAy:^yBy.A-l)('^ = J dW^'iy)dW^'iz) I dkj^^ [-lR\{yo)z'dl 

+ lRV,iyo)z'-%d^ - ci?(,o)) / d'if^ (30) 

Integrating by parts and removing the coordinates z through z idq, we can rewrite this expression as: 



iA-^\ByzA-y\)^'^ = J d^zdpdq 



— TV? \ {q ^ ^ m? 3 {{q — p)"^ — 
r(2-L>/2) /I 



(47r)^/2 (to2)2-D/2 



e Riyo) (31) 



In the last step we have done the z integration first and then those corresponding to the momenta. Finally integrating 
in m2 we find the effective lagrangian contribution to order 0{d^): 
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42(2/0) = + 1) ( J - R(yo) (32) 



Now we derive the divergences with four metric derivatives. In this case, there are contributions from A ^BA ^ 
and A'^BA~^BA-^: 



- 2^^Vft,,(2/o) - ^Raxiyo)R^\{yo) + ^K'f3xiyo)R\'',iyo) + ^Rj'f3,iyo)R\-,{yo) 
+ ^Kpniyo) + ^Ro.p-,iyo?l z'^z^z'di + (^li?^/^^,jyo) - YhR;Ay^)R\''M 

X z^z^z'z'^dldt - s}-K^p[y,)z'^zP ^fl^g-y^{z) I dq^^ (33) 
Removing the coordinates as before and integrating in positions and momenta we find: 

(Ay^yBy.A^y'o)''' = -J^^^^^ {-^QR''"'iyo)R,.M 

- ^R'-'iyo)R,Ayo) + J (J - ^) °^(yo)) (34) 

Notice that this resuh is finite. The divergences wiU appear when doing the ui^ integration. We add to this term the 
lowest order contribution from the next one, namely: 



^ / dW^^iy)d^zg'/'{z)dHg'^\t)dW^\u) / dkdqdp 



X (-lRMyo)z'di + lRV,iyo)z'z^d;d:;-iR{yo)^ ^yA.g^y^^,f.r^ 

X {-\R\{y.)u^dl + ^-Rr,{yoXu^d;d: - ei?(yo)) ^.9-^/^(«)^f^ (35) 
In a similar fashion to the previous case we find the results for the regularized integrals: 

iAy„\By.A-,^B,^AZy\r' = -J^^,^ {I e)'i?(j/o)^ - Y^i?-(2/o)i?..(yo)) (36) 
Substracting (^) from ( ^6| ) and integrating in we obtain the divergent lagrangian up to 0(9'*): 

>ci2(yo) - ^ (^■^R'"'"'iyo)R,.xpiyo) - ^R'"'iyo)RM - ^ (J - ^) °^(yo) + " ^^^°^') ^^^^ 

Comparing with the well-known Schwinger-DeWitt expansion, we see that C\^[j{yn) is proportional to ai(0,j/o) sjid 

'^dit)(yo) to 02(0, yo)- As far as the above expressions are scalars they will have the same form in any coordinate 
system, not neccessarily geodesic. We can then perform a coordinate change and integrate to obtain the corresponding 
EA: 

W[g^j^^]div^ J d*xCdiv{x) (38) 

We have included in Cdivix) the g^^'^{x) factor coming from the integration measure. It is possible, in principle, that 
when doing the coordinate change, new non-covariant terms appear provided they vanish only for geodesic coordinates. 
However, in absence of gravitational anomalies the EA is scalar and accordingly such terms are not permitted. On the 
other hand, terms with an odd number of derivatives yield terms with an odd number of momenta in the numerator 
which vanish in dimensional regularization. The above are the only possible divergences, higher derivative terms give 
rise to momentum integrals with more momenta in the denominator that turn out to be finite. 
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Apart from the divergences, which are purely local contributions, the EA also contains finite nonlocal pieces that are 
responsible for the pair creation processes. We have seen that in the QED case (||), the massless limit is well-defined 
due to the cancellation of the mass logarithmic dependence associated to the divergences with that coming from the 
nonlocal pieces. This fact allows us to extract some information about the nonlocal terms from the knowledge of the 
divergences. In next section we will profit this connection to find part of the nonlocal structure of the gravitational 
EA by means of a point-splitting procedure in the divergences. 



D. Nonlocal contributions 



In the previous section we have obtained local contributions up to order 0{Ti?). If we continue the calculation to 
higher orders, we would get a power series with terms of the following type: 

where V denotes the covariant derivative. This is a typical derivative (or adiabatic) expansion which is only valid at 
low energies. 

All the terms in (|3^ ) can be classified by the number of curvatures they contain. Those terms with a fixed number 
n of curvatures will give rise to a series with an increasing number of covariant derivatives. If we could add all these 
terms together, we would obtain nonlocal contributions that would provide the exact n-point Green functions (with 
curvatures in the external legs) [p^Pj9[|. These Green functions will be valid for any value of the mass m, as in the 
QED case (|). 

In this section we propose a method that effectively carries out that resummation for the quadratic operators, i.e 
we are interested in those terms of the form V"7?,^. The basic idea is to perform a point-splitting in the quadratic 
parts in curvatures in ( p3|) and (|35|), at the end of the section we will argue that this procedure gives rise to the correct 
resummation up to 0{TrPTZ^) terms. 

As we have just commented, all the terms in ( |39| ) are local and finite. The reason why we have not obtained nonlocal 
terms as in (||) is that the normal coordinates expansions are performed around a single point j/o- In (^) and ( |35| ) 
there are products of curvatures evaluated at the same point, i.e, TZ{yo)TZ{yQ) . Using again the normal coordinates 
expansion we can rewrite these products as: 

7^(2/o)7^(yo) = 7^(2/o)7^(y) + (7^(yo)V7^(yo)2/ + n{yo)W^TZ{yo)yy + ...) 

+ {n{yo)n{yo)'Jliyo)yy + niyo)Vn{yo)n{yo)yyy + ...) + ... (40) 

This expression allows us to split the points, the price to pay is the modification of the coefficients of the infinite 
higher order terms. We will not modify the linear part in TZ in (p^). Let us first obtain the nonlocal result and then 
we will argue that the new terms V"??,^ generated in (|4^) will exactly cancel those coming from the expansion in (|2^). 
By means of the above point-splitting, the TZ^ contributions in (^3|) and ( |35| ) will remain as: 

dq ^ -f 0{V^n^) (41) 



and: 

{A-^^By^A-lBtuAz')^-^'^ = / dW'\y)d^zg'/\z)dHg'/\t)dW'\u) I dkdqdp 



-lRMyo)z''di + \R^'',{yo)z^z^did: - myo)) '-^f^'^^f-^^ 

-\R\(t)u'>dl \r\ %{t)u^u^d-d: - mt)) ^p^f-''\u) + 0[V^n^) (42) 

where, as mentioned before, 0{\/^Ti?) denotes local finite terms with two curvatures and an arbitrary even number 
of derivatives. There are in principle different ways of performing the splitting, depending on the choice of the pair 
of points, but all them are equivalent up to higher order terms as can be seen from (40). 
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Removing the explicit coordinates occurrences, we obtain: 




Rivo) 



,2^,2 



3 [q^ ~ m 



2\3 



2 Rt^-{t){p + q)^ip + q ), 
3 ((p + (7)2-m2)2 



(43) 



Using the equations (pO|), (61) and ( p^ ) from the Appendix, neglecting higher order terms and integrating in m^, we 
obtain the following contributions to the effective lagrangian: 



1 



327r2 



dHdpe'P"^ {A- F{p^;m^)) 



1 /I 
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a Riyo)R{t) 



1 
108 



i?M^'(2/o)i?^''(t) 



o(v27^2 



(44) 



F{p^;'m?) being given in (||). 

Apart from the local divergent and nonlocal finite terms, in the dimensional regularization procedure, finite local 
terms do arise. However, their coefficients will be absorbed in the definition of the renormalized parameters and they 
will not be explicitly considered. 

In the same form as before we obtain from (p3): 



(yo) 



1 



327r2 



dhdpe'P' {A~ F{p^;m^)) 



0(V27^2) (45) 



Adding both contributions in ( [44[ ) and (|4^) and including the and 2 derivatives divergent contributions given in 
( p9| ) and (|3^), we can write the nonlocal EA in a slightly different notation: 



m 
~2 



- 771^ (A + 1) ( i - e ) R{x 



i-i?^''^''(a:)i?^.Ap(x) 



—R^-(x)F{U-vfi')R^,{x) , 2 



1 /I 



6 



-n R{x)F{U-^m^)R{x)\+0{m^Vf) 



(46) 



Here the F{^\vr?) operator action should be understood through the expressions (Q) y ( p5| ) as we did in the QED 
case (H), in addition we have neglected local finite terms and total derivatives. On the other hand, we see that the 
quadratic divergences agree with those obtained in the previous section. 

Consider now the massless limit of the EA in (46). For small masses compared with p, the Mandelstam function 
behaves as shown in (||). As in the QED case, the term proportional to log(m2), associated to the divergences in (|4^), 
has the same coefficient but with opposite sign as that coming from the Mandelstam function. Therefore they exactly 
cancel. This allows us to obtain the regular massless limit: 



180 



Rf^''ix)Tia)R^Ax) 



1 /I 



2 V6 



-n R{x)r{n)R{x) +oin^) 



(47) 



where r(n) is given in (||). This result is what one had expected on dimensional grounds however we have got in 
addition an explicit representation, completely analogous to the flat space-time case, for the form factors. In principle, 
the result in (^ cannot be considered, strictly speaking, as the 0(7?.^) contribution in the curvature expansion since 
in the calculation we have not considered those terms with two curvatures and an arbitrary number of covariant 
derivatives, denoted by 0{W^TZ^). We have obtained a nonlocal expression, but the infinite 0{W^TZ^) terms could in 
principle modify it. However, we have seen in the previous section that the divergences have a logarithmic dependence 
in (whose calculation is unambiguous). In addition, we showed in Section 2 in the QED case that the coefficients 
of the log(m2) in the divergent part have to be the same as those of the log(p2) in the nonlocal pieces, in order to 
have a regular masless limit. In ( |47| ) we have seen that our result indeed posseses a regular masless limit. Since the 
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Oiy^Ti?) terms are finite, they could only modify the finite local and nonlocal part, but not the divergences. However 
such a modification in the nonlocal part could spoil the regularity of the masless limit. Accordingly, we conclude that 
the Oiy^Ti?) terms can only modify the local finite pieces and those proportional to m? , but not the nonlocal ones 
containing log(m^). For that reason, the result in (^) includes all the quadratic curvature contributions except for 
those containing an arbitrary mass power denoted Oim^Tl'^). 

A compatible result has been obtained using partial resummation of the Schwinger-DeWitt expansion |^ and also 
by means of the so called covariant perturbation theory ||^, in this case it was possible to derive the cubic terms in 
curvatures for asymptotically flat manifolds. Notice however that the boundary conditions that normal coordinates 
impose on the metric tensor are slightly more general than the asymptotically Minkowskian condition. In fact, we 
have only required that the curvatures and their covariant derivative vanish at infinity. This condition includes, for 
instance, those Robertson- Walker manifolds in which the expansion rate a/a asymptotically vanishes in the future, 
although the scale factor itself a{t) does not tend to a constant. This kind of manifolds are known to accept the 
definition of the so-called adiabatic vacua In the possibility of defining such vacua was extremely useful for 
the interpretation of the effective action as vacuum persistence amplitude. 

To summarize, the procedure we have just presented makes it possible a partial resummation of the higher order 
terms in the EA, within the mentioned limits. Neglecting Oiji?) terms implies that ( ^6| ) will be a good approximation 
for VV7?. >> TZ^ . In a similar fashion, when we neglect 0{m?Ti?) terms we are assuming that WTZ » rn^TZ. 
However, if we are only interested in the two point Green functions with external curvature legs, the massless limit in 
(47) is exact. The renormalization of the EA can be done following the standard procedure in the classical references ( 
[ 13 and |Q ) . Here we will only mention that from the viewpoint of the Appelquist and Carrazone decoupling theorem 
[ 20 1 , the scalar field does not decouple from gravity since there are new terms in the EA ( ^6| ) which are not present 
in the Einstein-Hilbert action and they are not suppressed by powers of the particle mass m. 



IV. INTEGRATION OF THE STANDARD MODEL MATTER FIELDS: THE EA FOR TORSION 

Up to now we have only worked with scalar fields. However, in order to include the effect of the matter content 
present in the SM in the gravitational EA, we have to deal with the integration of fermionic fields. In this case we can 
proceed in a similar way using ( p^ in order to obtain the EA for the gauge fields and gravitation by integrating out 
the SM matter fields. The main novelty is that, in addition to the gravitational field (the vierbein), fermions couple 
also to the pseudotrace of the torsion as it was discussed at the beginning of the previous section. The SM matter 
lagrangian in a curved space-time with torsion can be written as [pi| : 

Cm = Va (Q(* - M'^)Q + C{t - Af^)£) (48) 

where: 

pQ ^ = r{d^ + + s^i,), p"^ - -f^Dj: = rid, + + ^^75) (49) 

with M*^ and the mass matrices of quarks and leptons. We have followed the notation in In the following 
we will concentrate only in the gravitational couplings, so that we have neglected the gauge fields contributions in the 
operators. 

In order to obtain the torsion contribution to the EA up to 0{S^), we first consider a flat space-time with torsion 
and afterwards we will include the effect of curvature. From the above lagrangian we see that torsion behaves as the 
electromagnetic field in QED, without considering the 75 coupling, which does not affect the calculation of terms with 
an even number of fields, in the massless limit. Using the result in (|^) and changing e ^ — > —1/8(^75, we obtain in 
this limit: 

W[S] = -(^^+Nnf + N,N,^ / rf'^^^'^r(n) V + 0{S^) (50) 

where r(n) is defined in (||), Nf is the number of families, N^, the number of neutrinos, N^f that of Dirac fermions 
and Nq the number of quark flavors. 5'^,^ = df^S^ ~ d^S,. Due to the absence of right neutrinos, there could be, in 
principle, parity violating terms. However it can be shown that those terms are total derivatives. 

In order to introduce the space-time curvature, we recall that the SM matter sector is locally conformally invariant 
(for massless fermions), this is also the case of the counterterms (see | 
has been studied for manifolds with torsion). When including the curvature, apart from those in (|50|), there could be 
quadratic terms in torsion in the generic form TZS^ . However such terms are not conformally invariant. Therefore the 



23 3] where the renormalization procedure 
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ones obtained above are the only possible divergences. Next we will confirm these results with an explicit calculation 
of the divergences in a space-time with curvature. 

The SM Dirac operators in Euclidean space are not Hermitian because of the electroweak gauge couplings and the 
absence of the right neutrinos. However, in Euclidean space the EA divergences are real [|l8| , thus it is enough to 
calculate: 2Re W[e, f ] = - log det(O), where O = {p + M^ip + M) and p and M denotes the joint Euchdean Dirac 
operator and mass matrix for quarks and leptons. The heat-kernel expansion together with dimensional regularization 
allows us to obtain: 

- iTrlog(O) = -^-^trf; [n - J d^x^g a^iO^x) (51) 

The well-known HMDS coefhcients a„ are given for the above operators by: 

a^iO.x) = 1, ar{0,x) = \r-X 

6 

a2(0, X) = ^[D^,D,][D^^, D''] + X]] + \^X^ - hiX 

where, when we neglect the gauge fields contributions, the operator can be written as: 

= D^D^ + X + M'', with X = 755';^ + 25^5^-i[7^7'^][d^,d,] 

and 

Writing the result in Lorentzian signature we have (see also the previous works IIJI): 



[e, f ] = - ^ ^ 7V„ / d^x 
( 1 7 



«.'l^ + |)-M?(A. + l)(-lfl + ls' 



1 „ „„., 1 



-|- A, f S^^ Sinj ^ RiiunrrR^^^^ ^ RiiijR^^ H R^ 

\ZM ^ 1440 ^ 180 ' 288 



(53) 



J2f denotes sum over the different families and J2i — J2iept ~^^cJ2quarks different flavors in each family, 

including quarks and leptons. On the other hand, Ni is the number of spinor components: 2 for neutrinos and 4 for 
the rest of fermions. Mi are the different fermion masses, = TV^ -I- \og{fi^ / Mf) and /i the renormalization scale. 

This formula is compatible with the flat space-time result in (|50|). There is no TZS^ term, as commented before, and 
we have discarded total derivatives. Following similar steps as for the scalar case, i.e. including log(n/Af^) factors 
and taking the massless limit, we obtain the flnite nonlocal contributions depending on curvature and torsion from 
the divergences. In the mentioned limit we have: 



1440 ^^-^^ ' ^ - 180 

where r(n) is given in (^) and 0(5*^) denotes terms with 3 or more torsion fields. Finally, due to the presence of 
chiral fermions, the EA could contain an abnormal parity sector, responsible for the gauge and gravitational anomalies. 
However, as far as the previous result has been derived from the real part of the EA (with normal parity), this sector 
is not taken into account. 

The renormalization procedure in this case will require, not only the introduction of quadratic terms and a constant, 
but also a kinetic and mass terms for the torsion field. Accordingly, the starting classical action should be: 

Sg^ J dx^ (^6^ + a^R'^-'^PR^^xp + a^R'^'^R^, + a^R^ + a^S^^S^^" + a^S^^ (55) 
The divergences can be absorbed by constants redefinition: 
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Uj a' aw(/u) — CjA^e + finite constants, j — 0..6 (56) 

where ao — ~2A/ {WttG) and ai ~ l/(167rG'). The values of Cj are shown in Table 3.1. 

The scale dependence of the renormalized constants is derived from the renormalization group equation [^,^ , in 
an analogous way to the scalar case. As a consequence, if at a given scale /i, the constants have certain values a^ip), 
their values at a different scale /x' will be given by: 

<(M') = <(M) + ^log(^) (57) 

where the Ci constants are those appearing in Table 3.1. In particular, the renormalized Newton constant G^{fi) has 
a scale dependence and therefore its value should be specified for a given fi. Thus G will depend on the size of the 
system we are considering and this could have an enormous importance in cosmology (see |25|1 ). 



Co 




C2 


Cs 










Ci 










48 / c 64^^/ 







Table 3.1: Renormalization constants for the SM case 

In the renormalized EA we can generate a torsion kinetic term by means of a finite normalization as follows: 
5^ — ^^'^{ii)S^, where a^ifJ-) = —■jZ^^{^j,). Thus in the massless limit we have: 

y^r^gr^ ^ J ^4^^ l^^^^'-p-' _^ nonlocal terms + 0{S^)^ (58) 

As a consequence the physical torsion field SJ^ will behave as an abelian gauge field. Therefore we have generated a 
kinetic term for torsion even starting from a theory without propagating torsion. 



V. CONCLUSIONS AND DISCUSSION 



In this work we have dealt with the computation of the low-energy Effective Action (EA) for gravity obtained when 
matter fields, both scalar and fcrmionic, are integrated out. As a much simpler exercise, we have started with by 
reviewing the low-energy EA for the electromagnetic field obtained when the electronic field is integrated out. Special 
attention has been paid to the nonlocal terms which are related to the particle production probabilities. 

In order to integrate the scalar fields to compute the EA for gravity, we have used the normal coordinate expansion. 
This method has been used previously to compute the divergent terms of the EA. By comparison of this kind 
of expansion around different space-time points, we have been able to extend the previous work to obtain also the 
nonlocal finite terms of the EA up to quadratic terms in the curvature. In particular we have arrived to the conclusion 
that the scalar field does not decouple from the EA in the large mass limit in the Appelquist-Carrazone sense. Finally 
we have also discussed the meaning of the expansion in the massless limit. 

In order to be able to consider the matter field content present in the Standard Model we have also studied the case 
of fcrmionic matter. The main novelty in this case is that this kind of fields can also couple to the torsion pscudotrace 
in addition to the gravitational field. Thus we have obtained the low-energy EA for the gravitational field and the 
torsion including the divergent and the nonlocal finite terms up to quadratic terms in the curvature and the torsion. 
We have discussed the renormalization of this EA and we have found that a kinetic term is generated for the torsion 
pseudotrace. This result is quite interesting since in the standard Einstein-Cartan action ||l|, the torsion does not 
propagate. 

Finally we would like to stress that the results obtained in this work can be useful for the study of some interesting 
physical effects. In particular, the low-energy EA can be applied for the study of the quantum stability of classical 
solutions of the Einstein equations of motion. More specifically, the nonlocal terms of the EA which have been 
computed here, are complex in general. Classical solutions that give rise to an imaginary part when the EA is 
evaluated on them, are unstable by particle radiation even if they are stable at the classical level. Moreover, the 
imaginary part of the EA evaluated on classical solutions, can be used to compute particle production rates and 
spectra thus providing an alternative method to the more commonly used based on the Bogolyubov transformations. 
The expression that we have obtained for the nonlocal form factors is particularly appropriate for the calculation 



12 



of particle production probabilities in cosmological space-times. The boundary conditions imposed by the normal 
coordinates expansion, allow to apply this method not only to asymptotically flat manifolds, but also to some other 
cosmological manifolds where the expansion rate asymptotically vanishes. Work has been done in this direction and 
it will be presented elsewhere fill . 
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VI. APPENDIX 



Dimensional regularization formulae 

In dimensional regularization |2^] , the space-time dimension is taken to be D = 4 — e and the poles are parametrized 
through Ac that was defined in Section II. We will use the notation dq = /.I'^d^q/ (27r)^ where fj, is the renormalization 
scale. Some useful expressions are: 

{q'Y ^ . (-1)'-"^ r(r + D/2)T{m ~ r ~ 
^[q^-R^]'" * (47r)^/2 r(D/2)r(m)(i?2)™-'-i)/2 ^ 

h{p,m')^ I dq—^ , ^^^tl)^r{3-D/2) C dt{l - t){^m^ + pH{l - t))^'^-^ 







((j2_m2)2((g+p)2_TO2) (47r)^/2 

( 

2(47r)2 • - ^: 



h [P, m^)dm^ ( + 2 - log ^ - Fip^; m') ] (60) 



^2 iP,m ) = / dq 



I2" {p^ m?')drn? 



(g2 — rn?Y{{q -\- p)2 — m?) 
*^(j^r(4 - D/2) j'^ dt{l tf{-m' + pH{l t)f'^-^ 



P^P"" + 9' 



^^-m^+pH{l-t) 



2(3- D/2) 



12(47r)2 



A, + ^ - log ^ - F{p';m') ] g^^ + 0{p') 



(61) 



J^..o.^3| 2. I ,~ q'^q-'q^q^ 
I3 ' {p,m ) = / dq 



((72 — m2)4((g + p)2 _ j^2^ 

^ ^' T{b-D/2) dt{l-tf{-m'+pH{\-t)f'^-^ ^ 9 +9 9 +9 9 



6(47r)^/2 ' Jo {6-D){8-D) 

+ 0{p') 

^3 iP^^)d-^ =96(4^r^+3"^°S^-(,'-^j^(^'™)-^j 

X (g"'^/'^ + .g"^g^'' -f g^^g^"") + 0{p^) (62) 

where F{p'^;M'^) is given in (|). These expressions have been obtained in the Minkowski space-time. We have not 
explicitly included the +ie terms accompanying the momenta in order to avoid the confusion with that coming from 
dimensional regularization. 



Normal coordinates expansions 



We will show the Riemann normal coordinates expansions for different objects up to order 0{d*). We will take the 
point yo as the coordinates origin. For the inverse metric tensor we have: 
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20 



(63) 



The metric determinant g = \ det g^i^j has the following expansion: 



g{y) = 1 + ii?o/5(yo)2/"/ - ^i?a/5ny"/y^ + ( ^R^pR^s - ^Rxcp.R^^s^ + :^i?o;3;75) 2/"/y'^2/' + ^(5^) (64) 
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For the square root of the metric determinant and its inverse we have: 



g^'^y) - 1 + ^i?a/3(j/o)y"/ - ^R^p-r^y'^y^y'' + [^^Ro.pR,8 - ^Rxc.fi.R\s^ + ^R^Pn^ y'^y^y^y' 

+ 0(9^), (65) 

5"'/'(?y) = 1 - ^i?aMj^0)y"/ + ^i?a/3;7y"/y^ + {^Ro.pR^6 + Y^i?Ac./3«i?\5'' - ^^a/3;7*) 

+ 0(9^) (66) 
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